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Superconducting heterostructures with spin-active materials have emerged as promising platforms 
for engineering topological superconductors featuring Majorana bound states at surfaces, edges and 
vortices. Here we present a method for evaluating, from a microscopic model, the band structure of a 
semiconductor film of finite thickness deposited on top of a conventional superconductor. Analytical 
expressions for the proximity induced gap openings are presented in terms of microscopic parameters 
and the proximity effect in presence of spin-orbit and exchange splitting is visualized in terms 
of Andreev reflection processes. An expression for the topological invariant, associated with the 
existence of Majorana bound states, is shown to depend only on parameters of the semiconductor 
film. The finite thickness of the film leads to resonant states in the film giving rise to a complex band 
structure with the topological phase alternating between trivial and non-trivial as the parameters 
of the film are tuned. 



A topological superconductor is a superconductor with 
a non-trivial topological order and associated zero en- 
ergy edge-, surface- or defect modes. Recently^ a clas- 
sification of the possible distinct topological phases of 
superconductors (and insulators), characterized by an 
integer topological invariant, was constructed following 
an earlier classification scheme based on the symmetric 
spaces of Bogoluibov-de Gennes (BdG) Hamiltonians due 
to Altland and Zirnbauer^. One of the most interest- 
ing properties of the more exotic classes of superconduc- 
tors is that the creation operators associated with the 
zero energy boundary modes are the same as their an- 
nihilation operators, leading to topologically degenerate 
ground states. Such modes, commonly referred to as Ma- 
jorana fermions, constitute essential ingredients in pro- 
posals of schemes for fault tolerant quantum computing 
known as topological quantum computing^. The proto- 
typical example of such a systerrP, the p x + ip y supercon- 
ductors, may be intrinsically realized in the 2D layered 
material Sr2Ru04. However, since time-reversal symme- 
try is spontaneously broken, the p x + ip y and p x — ip y 
phases are degenerate which generally leads to forma- 
tion of domains with both chiralities, thus complicating 
the edge-state structure. Furthermore, the layered struc- 
ture, as well as the low intrinsic gap at higher temper- 
atures implies that Majorana bound states at e.g. vor- 
tices are only protected by a very small minigap in the 
mK range. These considerations contributed to the gen- 
eral excitement following several proposals for engineered 
topological superconductors, where heterostructures in- 
volving conventional s-wave superconductors and various 
spin-active materials are used to produce artificial lower 
dimensional topological superconductors^"^. While the 
topological invariant of the "intrinsic" p x ± ip y super- 
conductors can be evaluated directly from the Bloch- 
functions (or equivalently from the Bloch-Hamiltonian) 1 
of the periodic structure, the "engineered" topological su- 
perconductors are highly non-homogenous in at least one 
spatial direction which complicates a direct evaluation of 
the topological invariant. The situation is simplified by 
phenomenologically modeling the system by adding an 



effective "pair-potential" to the uncoupled Hamiltonian 
of the non-superconducting material to account for the 
proximity effect which opens a gap due to the mixing of 
electrons and holes in the Andreev reflec tion process. In- 
deed such a model can be justifie d 11 * 12 ^ by an approach 
originally due to McMillan^, where the hybrid structure 
is described by a tunnel model, and elimination of the 
superconductor leads to an effective Hamiltonian of the 
system with an effective " pair-potential" . While provid- 
ing a very elegant qualitative derivation of the model, 
it should be noted that this is still a somewhat crude 
(depending on the detailed description of the tunneling 
terms) model of the proximity effect which obscures the 
role of Andreev reflection and may miss crucial aspects 
of the gap dependence on thickness, barrier height etc. 

In this article wc develop a framework to examine the 
spectral properties of proximity induced surface states in 
thin films on conventional superconductors. The film is 
assumed to be non-superconducting but generally spin- 
active. Aiming to retain as much generality as possi- 
ble, the thin film is initially treated as a "black box", 
characterized by a "reflection matrix" R, which acts 
as a matching condition for quasiparticle wave func- 
tions decaying into the superconductor (Andreev bound 
states/deGennes-Saint James states). Solvability of the 
matching conditions provide a spectral equation which al- 
lows us to obtain the effective 2D bandstructure e n (feii) 
(translational invariance in the directions parallel to the 
film is assumed) .The reflection matrix is a 4 x 4 unitary 
matrix with diagonal entries in electron-hole space. 

= H( fe||)£ )=(^ r J ib£) ). 

(1) 

The matrices are subject to symmetry constraints, in par- 
ticular the inherent particle-hole symmetry originating 
from a redundant description of the system 

r h (k\\,s) = r^(-k\\,-£), (2) 

as well as possible "true" symmetries such as time- 
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reversal symmetry, 

r e (fe||,e) =a y r^(-k\\,£)a y , (3) 

and spin-rotational symmetry 

r e (fc||,e) = r e (fc||,e)l, (4) 

with |r e | 2 = 1 due to unitarity. As will be shown in the 
next section, solvability of the matching condition of the 
exponentially decaying solutions on the superconducting 
side leads to an expression: 

= dct (1 - r- x a v e-^r h (T v e-^) , (5) 

with 7 = arccos(e/A) the Andreev reflection phase shift. 
The occurrence of the a y matrices reflects the nature of 
the Andreev reflection mechanism in a singlet supercon- 
ductor (e f— > h I and e J,— > —hf). The solutions of the 
spectral equation yields the 2D bandstructure s n (kn). 
We derive a representation of the spectral equation that 
allows for an intuitive interpretation of the spectrum in 
terms of quantization rules akin to the semiclassical Born 
quantization rules, ft should be emphasized, however, 
that this calculation is exact at this level. This represen- 
tation shall be used extensively in section [TT] where the 
reflection matrices are classified ac cordi ng to the scheme 
developed by AltlancP, and otheriPS]. There we shall 
also make general observations of the various symmetry 
classes and consider some representative examples which 
we analyze numerically as well as analytically to some 
detail. 

The central objective of this work is to find a micro- 
scopic description of the class of hetero-structures involv- 
ing spin-orbit coupled semiconductors with additional 
spin-splitting, originating either from an applied mag- 
netic field or proximity induced exchange coupling, and 
conventional superconductors believed to have topologi- 
cal^ non-trivial bandstructure for a certain regime of pa- 
rameters, and associated Majorana bound states at the 
edges of the systerrP^. For this class of systems, the 
aforementioned spectral equation reads: 

=cos 2 #sin(7 - 8x+) sm(j - 8X-) 

+ sin 2 6>sin(7 - £xV) sm (7~ <$X-)i 

where 7 = arccos(i?/A) is the phase shift from the An- 
dreev reflection process, 2Sx± — Xe± — Xh^ represent 
the differences in scattering phase shifts for electron-hole 
branches off opposite spin/hclicity (in this notation s 
is for singlet-like), while 25 X ± = Xe± — Xh± represent 
the differences in scattering phase shifts for electron-hole 
branches of the same helicity branch (t for equal-spin 
triplet-like). The rotation angle, 9, parametrizing the di- 
agonalizing matrix, U, describes the relative magnitude 
of spin-orbit coupling vs exchange splitting (6 = in the 
absence of spin-orbit coupling and 9 = n/2 in the absence 
of exchange splitting). 

The analysis of the spectral equation reveals that the 
system is always fully gapped for fen 7^ as long as spin- 
orbit interaction is present. At fen =0 (where spin-orbit 




FIG. 1. Schematic illustration of the system considered in 
this article: A bulk superconductor with a film of finite thick- 
ness, L, deposited on top. Translational invariance parallel to 
the film is assumed, while the Hamiltonian in the supercon- 
ducting region differers from that in the film region in that 
the superconducting order parameter A is nonzero only in 
the superconductor, while spin-dependent terms, here repre- 
sented by a term b(feii) ■ tr in the normal state Hamiltonian, 
are nonzero only in the film region. An interface potential, in 
the following modeled as a delta-function with weight Vb, is 
also assumed to be present. 



interaction is zero in any case), the gap opens and closes 
upon tuning the exchange splitting. The opening and 
closing of the gap signals a topological phase transition, 
in the sense that two distinct phases (bandstructures) 
can not be adiabatically (smoothly) connected without 
closing the gap (at critical point in parameter space cor- 
responding to the topological phase transition). The 
topological invariant characterizing the different phases 
is shown to be 

M = sign [cos(Sx°] e {-1, 1} ~ Z 2 , 

where 5 X ° = S x s + {e = 0,k { \ = 0) = -8xt(e = 0,fe|| = 0). 
Physically, a topological phase transition occurs when- 
ever one of the branches is pushed below zero energy 
which thus changes the number (more specifically the 
parity) of occupied bands. This interpretation is consis- 
tent with the recent work by Potter and Lee^, as well 
as Stanescu, DasSarma et. alP^ where they study the 
influence of multiple bands (due to quantized transverse 
modes) starting from the phenomenological model of the 
system. In our description multiple bands are an intrinsic 
feature of the microscopic model and the finite thickness 
of the film. 



I. SPECTRAL EQUATION 

The type of system in mind is a conventional (s-wave, 
singlet) superconductor with a non-superconducting film 
of arbitrary thickness deposited on top of it. Anticipat- 
ing the film to break spin-rotational symmetry we ex- 
tend the usual two-component Nambu spinor to its four- 
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component form 



4>i 



(6) 



where, by slight abuse of notation, we denote ip* — (ip T y . 
Due to the redundancy introduced in this notation, the 
four- component spinor is pseudo-real in the sense that 
= r x ty* . This structure of the spinor implies a con- 
straint on the Hamiltonian 



(7) 



usually referred to as particle-hole-symmetry. By this 
constraint the Hamiltonian then has the general structure 



H 



n a 

-A* -W 



(8) 



where H is the normal state Hamiltonian and A aa t — 
—g(4>a-ipa') represents the superconducting pair potential 
that should, in principle, be evaluated self-consistently. 
The attractive interaction —g is non-zero only in the su- 
perconducting region, implying that the pair-potential 
is zero in the film region, while the assumption of a 
singlet superconductor fixes the spin-structure of the 
pair-potential in the superconducting region A aa > — 
(ia y ) acr >A. Together with the translational invariance 
along the film, the Hamiltonian in the superconducting 
region takes the form 



ier j, A 

-£iiH#* 



(9) 



where £\\(-id z ) = ( lf n ] - n\\ within = [i — 2m - Thc 
elementary solutions of this Hamiltonian are plane waves, 
$(fc||, k)e tkz with the envelopes satisfying the eigenvalue 
equation 



#(fc N ,A) 



£n(k)u + AidyV = eu, (10) 



with a wave vector determined from the condition 



£ 2 =tf(fc) + A 2 



(11) 



Since the right hand side is a fourth order polynomial in 
k we generally have 4 different solutions. For e 2 < A 2 
we must have Re(£|i(fe)) = 0. By direct inspection 
we see that if k is a solution, then so are — k, +k* 
and — k* . The wave functions corresponding to these 
wave vectors are either exponentially growing or decay- 
ing into the bulk of the superconductor. For a semi- 
infinite geometry only the exponentially decaying solu- 
tions are physical, and we must have Imfc < 0. Denot- 
ing the solution with positive real part simply by k + we 
note that fc_ = — k*L is also a physical solution, while 



the other two are non-physical (not normalizable) . Let 
us also introduce the phase 7 = arccos(e/A) such that 
COS7 = e/A and sin 7 = |£ N (fc)|/A = y/l - (e/A) 2 . 
With this parametrization Eq. (flOl) obtains the form 



(12) 



where ti* = u(k±) and v^ 1 = v(k±). The full wave- 
function on the superconducting side is a linear combi- 
nation 



(13) 



where v = Ke(k±)/m and we have associated the wave 
function with positive real part of the wave vector as an 
"incoming" wave and the wave function with negative 
real part of the wave vector as "outgoing" wave. Specif- 
ically we have 



in ,.' 



<|jOUt _ 



(14) 

Applying the boundary condition <I> out = _R<I> m between 
the incoming and outgoing wave functions through thc 
scattering matrix R from Eq. ([!]) we get 



u = r e u^ , 
(ia y )e' l ' y u~ = rh{icr y )e~ l ~'u + . 



(15) 



Combining the two finally leads to the equation 

= (1 - r- 1 a y e- t -<r h a y e-^)u + . (16) 

The solvability of this equation then implies the spectral 
equation 



- det(l - r- 1 {k ll ,e)a y e^'<^r h {k lb e)a y 



-17(e) 



(17) 

The solutions of this equation constitute the bound state 
spectrum, or 2D bandstructure e n (kii), of the states lo- 
calized around the film region. 

It is useful to point out that a gapped system can be 
identified by the condition (note 7(0) = 7r/2): 



O^det 1+r" 1 (fe, \,0)a y r h (k\ \,0)a y ), Vfe N . (18) 



II. SYMMETRY CLASSIFICATION 

In this section we shall discuss the properties of the 
reflection matrix r e . Note that by virtue of the inherent 
redundancy, here referred to as particle-hole symmetry, 
the properties of r e directly determine those of rv In 
order to systematically analyze the properties we follow 
the classification scheme of Altland and Zirnbauer. 
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Classification of Superconductors 


Class 


Time-rev. 


Spin- rot. 


D 


X 


X 


C 


X 




Dili 




X 


CI 




✓ 



Let us first discuss the effect of time-reversal symmetry, 
Eq. ([3]) which, together with particle-hole symmetry, Eq. 
([2]) implies a y rh(k\\,e)a y = r e (fc||, — e) and thus the con- 
dition ( fl8| for a gapless system reduces to 

det(l + r-- 1 (fc M ,0)T' e (fc||,0)) ^0 



VA; 



II' 



(19) 



which is always fullfillcd. Breaking of time-reversal sym- 
metry is thus crucial in order to observe any topolog- 
ical phase transitions (which occur at gap-closings) in 
proximity-based superconducting 2D systems. Since no 
topological phase transitions can occur, we note that the 
system must always be in the trivial phase. We can thus 
conclude that, for proximity structures, the only accessi- 
ble topological phase in symmetry classes CI and Dili is 
the trivial phase. 

In the context of Majorana bound states it is further- 
more crucial to also break spin-rotational symmetry since 
otherwise any quasiparticles in the system will have the 
form 7 = uifjl + vt/j^ or 7 = uij}\ + vip^ neither of which 
can possess the defining property of Majorana fermions 
7^ = 7. Thus it should be clear that the most interest- 
ing class (at least in the context of Majorana fermions) 
is the class D. Nevertheless, for pedagogical reasons we 
shall first consider the two classes CI and Dili in the 
next sections, followed by an analysis of the most inter- 
esting class D. While a complete analysis, including class 
C, would be desirable, it falls outside the scope of the 
present article. For interested readers we refer to the lit- 
erature on proximity structures with orbital coupling to 
the magnetic fields^"^. 



A. Class CI 

If the system preserves spin-rotational symmetry, the 
reflection matrix satisfies equation Q. Particle hole 
symmetry then requires = e %Xh \ with Xh(k\\,e) = 
% e ( — fen , — e). The spectral equation can then be written 

as 



(1 



-i(2 7 ( £ )+ Xc - X , l ) 



) = o. 



(20) 



Time-reversal symmetry requires Eq. ([3|, which implies 
Xe(k\\,e) = Xe(-fe||,e), andxh(k\\,e) = Xe(k\\,-e). The 
spectral equation then reduces to the the "quantization 
condition" 



2j(s) + Xe(k\\,e) - Xe(k\\, 



2im. 



(21) 



The physical interpretation of this equation can be un- 
derstood from Fig. [2] (Top). To obtain a qualitative 



understanding of the solutions to this equation we may 
employ the Friedel sum rule: 

Xe(k\\,e) - Xe(k\\, -e) = n J de'p(k\\,e'), (22) 

where p(k\\, e) is the normal density of states in the film 
on top of the superconductor. If we assume that there 
exists a single (spin-degenerate) resonant band £"o (fei 1 ) 
with a broadening T, i.e. the density of states has the 
generic form 



P{k\\,e) 



r 



(23) 



7rr 2 + ( £ - £o (fc||)) 2 ' 
then the phase difference has the form 
Xe(k\\,e) -Xe(k\\,-e) = 

1 arctanf £ + £ ° (fc '^ ^ 



Taking the cosine of Eq. ( 21 ) and expanding for small F 



one arrives at a BCS-like dispersion equation 
£ 2 =4( k \\) 



cfl'< 



A 



cir 



(25) 



The solution of the full dispersion equation is depicted in 
Fig. ([2]) for the case of a resonance-band with parabolic 
dispersion eo(k\\) = fcjj/2m — \i. The main feature, of 
course, is the hybridization between the electron- and 
hole-branches of the resonant bands due to the Andrccv 
reflection. The induced gap is given by A c ff = \/AT. If 
there are multiple resonant bands in the film the problem 
can be treated analogously provided the energy-spacing 
between the bands is larger than the width T. The An- 
dreev reflection will then induce hybridization between 
the electron- and hole branches of all bands (not neces- 
sarily at the chemical potential). 

While this is a convenient and intuitive description of 
the problem, it does not provide any information about 
the bandstructure in terms of microscopic parameters. 
To this end, we model the semiconductor film region by 
an effective mass Hamiltonian and the interface between 
the superconductor and the film by a (S-function potential 
of strength Vq. To account for the finite thickness L of 
the film we impose hard-wall boundary conditions on the 
other boundary of the film. The scattering phase shift is 
then given by Xe(k\\,e) = n + 25 e (k\\, e) with 



tan<5 P = 



v sin(fc'L) 



v 1 cos(fc'L) + 2y sin(fc'L)' 



(26) 



where v(k\\,e) = Re(fe(fe||, e))/m, v'(k\\,e) — 
Re(A;'(fe||, e))/m' denote the velocities and m, m' the 
masses in the superconductor and semiconductor regions 
respectively. The resonance condition k' L — nir cor- 
responds to resonance bands with energies e ni o(fe||) = 

2m 7 VTr) +2™"^ Naively, one might try to expand the 
left hand side of Eq. ( 26 1 around one of the resonances 
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FIG. 2. (Top) Schematic illustration of the Andreev reflection 
process: Transmission through the semiconductor film and re- 
flection at its surface produces an overall phase-shift Xe/h for 
electrons/holes. Resonances in the film are related to these 
phase shifts through the Friedel sum rule ( 22 1. Andreev reflec- 
tion on the superconducting side connects the electron- and 
hole branches through an additional phase 7. A bound state 
is formed through the cyclic process of electron reflection on 
semiconductor film - Andreev reflection - hole reflection on 
semiconductor film - Andreev reflection - etc. The dispersion 
equation ( |21[ | can be interpreted as a type of Born quan- 
tization rule. (Bottom:) Solution to the spectral equation 
|25| indicated by thick (red) line for weak (Left, V = 0.03A) 
and strong coupling (Right, V = 0.2A). The resonant bands 
are indicated by the parabolic thin and dashed curves (blue), 
while the effective gap A e fr is indicated by the horizontal 
dashed lines. Notice also the bending of the curves as the 
band approaches the bulk gap A. 



£n,o(fe||) to obtain a form similar to Eq. (|24t), and thus 
obtain a spectral equation like Eq. |25| ) . Inthis case the 
width of the resonance would be r(fe N ) = vL ^Mp 1 ^ , 
However, it turns out that the thusly evaluated effective 

gap, A c ff = yj^f-^rj^, neglects the effect of the inter- 
face potential oc Vo, and thus overestimates the gap. By 
instead considering the low energy limit of the spectral 
equation, by keeping terms up to second order in e, one 
obtains a much more realistic expression for the induced 
gap which in the limit 2mVoL max[fcp.L, ^ 2 mV l \ 
reduces to 



Act? 
A 



(27) 



With this expression for the effective gap, the expression 



e±(fe||) = ±Jel(k\\) 



crt ' 



(28) 



fits reasonably well with the numerical solution for the 
lowest energy band around k\\ = (see Fig. |3j) in the 




FIG. 3. 2D Bandstructure of a typical film of thickness 
L — 3A_f in class CI, i.e. respecting time-reversal and spin- 
rotational symmetry, are depicted with thick full lines, while 
the normal resonance-bands ±e n ,a(kit) are depicted by thin 
dashed lines. The thick dashed (red) lines show the disper- 
sion from Eq. 1 28 1 with the approximated effective gap (27 1. 



The fit works reasonably well for a weakly coupled system 
(left, Vo = vf) while for a strongly coupled system (right, 
Vo = 0.05uf) the fit becomes progressively worse eventually 
approaching the solutions ( 30 1 obtained in the strong coupling 



limit of the spectral equation, here indicated in the right fig- 
ure by the dashed (blue) horizontal lines. 



weak coupling limit 2Vq > vp. In the opposite limit, 
with v = v' , we can expand the spectral equation in A//z 
(Andreev approximation) to recover the deGennes-Saint 
James states^ satisfying the spectral equation 



— =±C»s (2tt^ V 



(Hi)' 



(29) 



with the approximate solution for the lowest energy 
states given by 



-1 



±A 



1 + 2(2^)2 



(30) 



where £0 = nvp/A is the superconducting coherence 
length. 



B. Class Dili 

For a system which respects time-reversal symmetry, 
but not spin-rotational symmetry, the reflection matrix 
is constrained according to Eq. (3]) which together with 
particle hole symmetry reduces the spectral equation to 



= detfl - r c 7 1 (fe||,e)e- i7(e V e (fe||,-e)e 



-17(e) 



(31) 



Let us assume, as is the case when the Hamiltonian in 
the semiconductor film contains a spin-orbit term, that 
the unitary matrix U which diagonalizes the reflection 
matrix r e does not depend on energy: 



r e (fe||, £ ) = C/(fe||)e i ^ fe H' e )t/t(fe||), 



(32) 
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where x = di a g(x+?X-) is a diagonal matrix contain- 
ing the scattering phase shifts in the diagonalized basis, 
i.e. the corresponding phase-factors are the eigenvalues 
of the reflection matrix. If this is the case, then the 
same unitary matrix also diagonalizes the reflection ma- 
trix r e (fe|i,— e), and the spectral equation reduces to 

= dct(l - e - l x(fe|h£) e -«7(£) e ^(fc||,-e) e - l7 (e)y (33) 

This equation is similar to the one in the class CI, with 
the important distinction that it decouples into two in- 
dependent spectral equations 



27(e) + X+(k\\,e) - x + (fe||,-e) = 2im, 
27(e) + X-i k \\,e) - x-(k\\,-e) = 2im. 



(34) 



This situation applies to semiconductor films with spin- 
orbit interaction, or more generally to any Hamiltonian 
of the form H(z > 0) = (-id z ) 2 /2m' — fi + b(k\\) ■ cr 
with 6(— feii) = — 6(fe||). Such a Hamiltonian can be di- 
agonalized by a unitary transformation which turns the 
vector 6(fe||) into the direction of the spin quantization 
axis. This unitary matrix then also diagonalizes the re- 
flection matrix, assuming the interface potential is spin- 
rotationally invariant. For each helicity quantum number 
a = ± we have separate sets of normal resonance bands 
in the semiconductor film given by 

^o(fe||) = ^ (f ) 2 + H + .|6(fe |l )|-M. 05) 

For the specific case of Rashba spin-orbit interaction due 
to the structure inversion asymmetry in the z-direction 
we have &(feii) = af- (k y , — k x , 0), while for Dresselhaus 
spin-orbit interaction due to bulk inversion asymmetry 
in the same direction (the growth direction) we have 

The spectrum of the 2D semiconductor can now be 
evaluated by following the same line of arguments as 
in section [II A| with the important difference that there 
are now two independent branches of opposite helicity 
a = ±. Since the spectral equation ( [34] ) does not mix 
states of different helicity, only electron-hole branches 
within the same helicity-branch will hybridize. The ab- 
sence of anti-crossing between electron-hole branches of 
opposite helicity is a direct consequence of the time rever- 
sal symmetry. Furthermore, two branches e n! +.o(fc||) and 
£n.-.o(fc||) cross at the time- reversal invariant momentum 
fen = as a consequence of b(— fen) = — 6(fc||) as re- 
quired by time-reversal invariance. In the weak coupling 
regime the lowest energy bands can, for small fe|| have 

the approximate form e±, CT (fe||) = ±^e^ tTfi (k\\) + A^ ff ) 
for the last normal state resonance band n which crosses 



the Fermi-level and the effective gap given by Eq. (271 



This spectrum is consistent with a BCS model of the form 
H = E*,,,* [^A k \\)4(k\\)cAk\0+AescUk\\)4(-k\\)] 
and can thus be considered as induced equal-helicity pair- 
ing. In the strong coupling regime the lowest energy 
bands again approach Eq. ( 30 1 (see Fig, [I| . 
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FIG. 4. 2D Bandstructure (thick black lines) of a semicon- 
ductor film of thickness L = 3Xf with Rashba spin-orbit in- 
teraction a^ Q = 0.3 A/ kp in the weak (left, Vq = Vf) and 
strong (right, Vo = O.Ivf) coupling regime. The normal he- 
lical resonance-bands are depicted by thin dashed lines (red 
and blue). 



C. Class D 

Before proceeding to the main objective of this article 
- a semiconductor film with spin-orbit and exchange in- 
teraction - let us first consider the subclass of D where 
spin-rotational symmetry is not completely broken. This 
is the case if we have exchange interaction in the semicon- 
ductor film, and the normal state Hamiltonian in this re- 
gion obtains a term oc cr z b z , where b z denotes the strength 
of the exchange interaction. Since spin-rotational invari- 
ance is preserved around the z-axis the reflection ma- 
trix is diagonal and takes the form r e (fe|i,e) = e ix ^ fc H ,e ^ 
with x = diag(x^, xi)- Due to particle-hole symmetry we 
have rh(ku,e) = e %x (~ h u~ e \ and the spectral equation 
becomes 

= dct(l - e-^^^crye-^e^-^'-^aye-^y 

(36) 

which again decomposes into two sets of spectral equa- 
tions 



27(e) + Xt( fc lh £ ) - Xi(-k lb -e) = 2mr, 
27(e) +Xl(k\\,e) - xt(-k\\,-e) = 2mr. 



(37) 



The effect of this spectral equation is that the electron- 
branch of the normal state resonance bands with spin f 
will hybridize with the hole-branch of the normal state 
resonance bands with spin J, and vice versa. Spin- 
rotational invariance around the z-axis implies that spin 
is a good quantum number. Here we should consider that 
a hole with spin-down can be interpreted as a particle 
with spin-up and vice versa. This leads to the formation 
of bands with spin-up and spin-down without avoided 
level-crossings between the two. If we write the normal 
state resonance bands as £ n ,cr,o(fe) = E«,o(fe||) + cb z , one 
finds that in the weak coupling regime the lowest energy 
bands have the approximate form 



£±A k \\) = ± y £ o( k w 



(38) 
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From this equation it is clear that the main effect of 
the exchange interactio n is to shift the energy levels 
up/down. At \b z \ = £q{Q) + A 2 ff » fi Q we reach a 
point where two branches of opposite spin cross each 
other at (feiue) = (0,0). This feature can be investi- 
gated more carefully from the spectral equation ( |37[ ) by 
setting e, fen =0 for which we get x CT (0, 0) — X-<x(070) = 
(2n — 1)tt. Applying again the Friedel sum rule we can 
write this condition as 



N t -N i = (2n - 1)tt, 



(39) 



where N a is the number of normal state electronic states 
at fen with spin a below e = 0. Thus the gap closing 
at this point coincides with the point where, in the nor- 
mal state, the number of electronic spin-up bands minus 
the number of spin-down bands below the gap changes 
modulo 2. This identification will be useful to us in 
the following when we combine spin-orbit interaction to- 
gether with an exchange field. To this end, we write the 
normal state Hamiltonian in the semiconductor film as 



-id* 



b(k 



IN 



M|| + b( k \\) ■ v, 



(40) 



The last term clearly breaks spin-rotational symmetry, 
but also breaks time-reversal symmetry since 6(— fen) 7^ 
— 6(fe||). The Hamiltonian can be diagonalized through 
a unitary matrix 



v = [cos §e-;v 



sm ■ 



cos 



(41) 



where (sin^cos^jSin^sintpjCOsfl) = 6(feii)/|6(fe|i)|. The 
same matrix also diagonalizes the reflection matrix 
r e (fe|i,e) and thus we can write r e = Ue tx U^ where 
again \ — diag(x+,X-) and Xa = ^ + 25 CT are defined 
analogously to ([26]) but with k' — > k' a for each helicity 
branch a = ±. Since time-reversal symmetry is broken, 
OyThUy can not be diagonalized by the same matrix U. 
From the only remaining symmetry - particle-hole sym- 
metry - we have 



r h (fc||,e) = f/*(-fe||) e ^(- fe H'- £ )t/ T (-fe||). 



(42) 



Using this parametrization and evaluating the determi- 
nant directly one obtains the following spectral equation 



= cos 2 #sin(7 - Sx s +) sin(7 - SxL) 
+ sin 2 6>sin(7 — 8x+) sin(7 — 5x-), 

where we have defined 



$x s A k \\i s ) 



Xa(fe||,e) 



Xa(-fe|h 



-e). 



(43) 



(44) 



The subscripts s and t indicate hybridization between 
electron-hole branches of opposite and equal helicitics 
and thus, in a sense, induced singlet-, or triplet super- 
conductivity respectively. Setting e = in this spectral 




^[TT]::::^: 



FIG. 5. Andreev reflection processes in the presence of spin- 
orbit interaction and Zeeman/exchange splitting. The normal 
state resonant bands are associated with the phase-shifts cor- 
responding to the eigenvalues of the reflection matrices r e 
and Vh- Due to the absence of both spin-rotational symme- 
try and time-reversal symmetry, diagonalizing r e , here in- 
dicated by straight lines in the intermediate (white) region, 
does not simultaneously diagonalize rh, thus depicted in the 
figure by crossed- and straight lines. Thus closed cycles in- 
volve electron- and hole branches, which become mixed by 
Andreev reflection, of both helicities. The mixing is deter- 
mined by the angle 8 in Eq. | |43[ ). Consequently there is a 
hybridization of electron-hole branches between all resonant 
states. The exception happens at feii,£ = where spin-orbit 
interaction is absent. 



equation we can identify points where the gap closes. Us- 
ing 7(e = 0) = ir/2 along with the fact that the scatter- 
ing phase shifts only depend on |b(fe||)| and everything 
else is isotropic, i.e. x CT (— fc||,e) — x CT (fc||,e) and thus 
8xl{e = 0) = and 8 x %(e = 0) = -5 X -^ = 0), we get 



: 0(fc||)cos 2 (x+(fc||,O)- X -(fc||,O)) +: 



: 0(fc||). 
(45) 



Since both terms on the right hand side are positive, the 
spectral equation can only be satisfied if both are equal 
to zero at the same time. This then reduces to the two 
constraints which have to be simultanously fulfilled: 

sin0(fc||) = O, cos( X+ (fc||,0)- X -(fc||,0)). (46) 

The first one is only zero at fe|| = we find that this 
is the only place where a gap-closing may occur. Fur- 
thermore, the second condition is the same condition as 
discussed in the beginning of this section and is thus de- 
termined by Eq. (391, i.e. the gap closes whenever the 



parity of the difference of occupied normal state spin- 
up/spin-down states is changed. Thus the gap closing is 
associated with a band-inversion and the different gapped 
phases can be associated with the integer number neZ 
from Eq. ( 39 1 . Distinguishing in particular odd, and 
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FIG. 6. (Left:) 2D Bandstmcture (thick black lines) of a 
semiconductor film with Rashba spin-orbit interaction = 
0.5A/A;f and exchange/ Zeeman splitting b z — 0.3A of thick- 
ness L = 3Af in the weak coupling regime. The normal helical 
resonance-bands are depicted by thin dashed lines (red and 
blue). (Right:) Gapped phases characterized by the Z2 in- 
variant M from Eq. | |47[ ) as a function of exchange splitting 
b z and thickness of the film L in the weak coupling (Vb = vf) 
regime. Here, beige refers to M. = while purple refers to 
M = l. 







k u /k F 



FIG. 7. 2D Bandstmcture (black lines) of a semiconductor 
film with Rashba spin-orbit interaction — 0.5A/kF and 
exchange/ Zeeman splitting b z — 0.1 A of thickness L = 3Af 
in the weak (left, Vo = «f) and strong (right, Vb = O.Ivf) 
coupling regime and the dispersion obtained from the effective 
model Eq. (48 1 depicted by dashed red lines. 



even parity we can define a Z2 number 
M = (N + -N_)mod 2 



3ign(cos( x+ (0,0)-x-(0,0))) G Z 2 . 



(47) 



This number, although not derived directly by topologi- 
cal arguments, can serve as a topological invariant in the 
sense that the number does not change within a gapped 
phase, only when the gap is closed can we go from odd 
or even parity. The corresponding bandstructure and 
"topological phase-diagram" is depicted in Figure [6] In- 
terestingly this number does not depend on the proper- 
ties of the bulk superconductor but only on the properties 
of the semiconductor film (of course a coupling to the su- 
perconductor must exist or else the formalism becomes 
meaningless). For weak coupling, the four lowest energy 
bands approximately follow the dispersion relation 



A c 2 ff 



f 2 



*so&||| 5 



±2 



(48) 



derived^ from the phenomenological model where a fi- 
nite pair potential is introduced by hand, with the in- 
duced pair potential given by the approximate expression 
( p7| ). The fit works well for weak coupling, but breaks 
down for stronger coupling (see Fig. [7J). 



III. CONCLUSIONS 

We have presented a general method of extracting the 
2D bandstructure of films deposited on top of a conven- 
tional superconductor. Using the symmetry classifica- 
tion of Altland and Zirnbauer we concluded, based on 
general arguments, that for classes CI (spin-rotational 
symmetry and time-reversal symmetry) and C (spin- 
rotational symmetry only) the only accessible topolog- 
ical phase is the trivial one, while in class D there exists 
a Z 2 number which can be expressed through the dif- 
ference of spin-up vs spin-down scattering phase shifts 
at the particle- hole invariant point fen = = e as 
M = sign(cos(x+(0, 0) — X-(0, 0))). The change in this 
number is accompanied by a closing of the gap in the 
2D bandstructure, in accordance with the typical notion 
of a topological phase transition. We furthermore de- 
rived analytical and numerical results for the 2D band- 
structure in materials with spin-orbit coupling and/or 
Zeeman/exchange-splitting. 
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